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PREFACE 


This  report  describes  a proposed  technique  for  predicting  overtopping 
rates  for  structures  exposed  to  irregular  wind-generated  waves  by  extend- 
ing the  method  of  predicting  overtopping  for  waves  of  constant  height  and 
period  presented  in  the  Shore  Protection  Manual  (SPM)  (U.S.  Army,  Corps 
of  Engineers,  Coastal  Engineering  Research  Center,  1975).  Results  pro- 
duced by  the  technique  should  be  more  representative  of  conditions  in 
nature  than  the  SPM  procedures.  However,  while  the  present  SPM  proce- 
dure is  based  on  laboratory  tests,  the  proposed  technique  has  not  been 
directly  tested  in  the  laboratory  or  the  field.  Therefore,  results  from 
the  technique  should  be  used  for  general  guidance  until  results  of 
laboratory  tests  presently  underway  become  available. 

This  report  was  prepared  by  John  Ahrens,  Oceanographer,  under  the 
general  supervision  of  Dr.  R.M.  Sorensen,  Chief,  Coastal  Structures 
Branch. 

The  author  acknowledges  the  extensive  computer  programing  assistance 
from  M.  Titus,  and  the  useful  suggestions  by  Dr.  R.M.  Sorensen. 

Comments  on  this  publication  are  invited. 


Approved  for  publication  in  accordance  with  Public  Law  166,  79th 
Congress,  approved  31  July  1945,  as  supplemented  by  Public  Law  172, 
88th  Congress,  approved  7 November  1963. 
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CONVHRSION  FACTORS,  U.S.  CUSTOMARY  TO  MJiTRIC  (SI) 
UNITS  OF'  MFASUIUiMCNT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  irches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  foot 

0.0929 

square  meters 

cul)ic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

mi  les 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

Knots 

1.8532 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

mi  1 libars 

1.0197  X 10" 

^ kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

ki lograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0. 1745 

radians 

Fahrenheit  degrees 

5/0 

Celsius  degrees  or  Kelvins^ 

*To  obtain  Celsius 

(C)  temperature 

readings  from  Falireniieit  (F)  readings. 

use  formula:  C » 

(5/9)  (F  -32). 

To  obtain  Kelvin 

(K)  readings,  use 

formula:  K « (5/9)  (F  -32)  ♦ 273.15. 

5 


SYMBOLS  AND  DEFINITIONS 


dg  water  depth  at  toe  of  structure  (feet) 

g acceleration  of  gravity  (feet  per  second  squared) 

equivalent  deepwater  wave  height  (feet) 

(iy)g  equivalent  deepwater  significant  wave  height  (feet) 

h height  of  the  structure  crest  above  the  bottom  (feet) 

ratio  of  Rp  to  Rg  for  a particular  probability  of  exceedance 

p probability  of  exceedance 

overtopping  rate  associated  with  a particular  probability  of 
exceedance  (cubic  feet  per  second  per  foot) 

%,  overtopping  rate  for  monochromatic  waves  (cubic  feet  per  second 
per  foot) 

a,  0^  empirical,  dimensionless  wave  overtopping  parameters 

Qj,  overtopping  rate  for  irregular  waves  (cubic  feet  per  second 
per  foot) 

Qq  j peak  overtopping  rate  (cubic  feet  per  second  per  foot) 

R runup  for  monochromatic  waves  (feet) 

Rp  wave  runup  associated  with  a particular  probability  of 
exceedance  (feet) 

Rg  wave  runup  with  deepwater  significant  height  and  period  (feet) 


6 


PREDICTION  OF  IRREGULAR  WAVE  OVERTOPPING 

by 

John  Ahrens 
I . INTRODUCTION 

This  report  gives  a proposed  technique  for  predicting  the  average 
and  peak  overtopping  rates  caused  hy  wind-generated  waves  acting  on 
coastal  structures.  No  guidance  for  predicting  the  overtopping  rates 
for  irregular  wave  conditions  is  currently  available  in  the  Shore  Pro- 
tection Manual  (SPM)  (U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering 
Research  Center,  1975)^,  or  other  Corps  publications.  Therefore,  this 
interim  guidance  is  provided  until  the  results  of  the  Coastal  Engineering 
Research  Center  (CERC)  laboratory  study  of  runup  and  overtopping  by 
irregular  waves  are  available. 

The  approach  is  consistent  with,  relies  on,  and  extends  the  SPM 
section  on  wave  overtopping.  In  addition  to  the  proposed  technique, 
the  method  provides  a simple  but  plausible  means  of  judging  the  con- 
servatism or  nonconservatism  of  the  traditional  approach  to  predicting 
overtopping  based  on  monochromatic  wave  conditions.  In  applying  the 
procedures  described,  a word  of  caution  is  noted:  Some  larger  waves 
in  the  spectrum  may  be  depth- limited  and  break  seaward  of  the  structure. 
In  this  case,  the  rate  of  overtopping  may  be  overestimated. 

II.  MONOCHROMATIC  WAVE  OVERTOPPING  EQUATION 

Equation  7.6  in  the  SPM  (Sec.  7.22)  is  given  to  compute  wave  over- 
topping. Weggel  (1976)^  provides  the  background  information  on  how  this 
equation  was  developed  and  gives  additional  information  on  its  use  in 
supplementing  the  discussion  in  SPM.  The  overtopping  equation  from  the 
SPM  is  written  as 


'0.217  , /h-dg\ 

L “ v " / 


in  which 


0 < 


h - d. 


< 1.0  , 


(1) 


^U.S.  ARMY,  CORPS  OF  ENGINEERS,  COASTAL  ENGINEERING  RESEARCH  CENTER, 
Shore  Protection  Manual,  2d  ed.  , Vols.  I,  II,  and  III,  Stock  No. 
008-022-00077-1,  U.S.  Government  Printing  Office,  Washington,  D.C. , 
1975,  1,160  pp. 

^WEGGEL,  J.R.,  "A  Wave  Overtopping  Equation,"  Proceedings  of  the  ISth 
Conference  on  Coastal  Engineering , 1976. 
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where 


g 

Hi 

h 

ds 

R 


“.05 


overtopping  rate  (volume  per  unit  time)  per  unit  structure 
length 

gravitation  acceleration 

equivalent  deepwater  wave  height 

height  of  the  structure  crest  above  the  bottom 

water  depth  at  the  structure  toe 

runup  on  the  structure  that  would  occur  if  the  structure 
were  high  enough  to  prevent  overtopping  (corrected  for  scale 
effects,  SPM  Sec.  7.21) 

empirically  determined  coefficients  that  depend  on  incident 
wave  characteristics  and  structure  geometry 


The  SPM  (Sec.  7.22)  gives  values  of  a and  for  a wide  range  of 

structure  types  and  wave  conditions. 

Equation  (1)  provides  an  unusually  good  fit  to  observed  overtopping 
rates  from  laboratory  tests  using  monochromatic  waves  (waves  of  constant 
height  and  period).  was  used  in  equation  (1)  rather  than  Q to 

specify  overtopping  by  monochromatic  waves.  The  following  section  in- 
cludes some  modifications  to  equation  (1)  to  adapt  it  to  irregular  wave 
conditions  so  that  overtopping  rates  for  irregular  waves  can  be  calculated 
and  compared  to  overtopping  rates  for  monochromatic  waves. 

111.  ADAPTION  OF  MONOCHROMATIC  WAVE  OVERTOPPING  EQUATION 
TO  IRREGULAR  WAVE  CONDITIONS 


The  fundamental  assumption  in  adapting  equation  (1)  to  irregular 
wave  conditions  is  that  the  runup  caused  by  these  conditions  has  a 
Rayleigh  distribution  of  the  type  commonly  associated  with  wave  heights. 
This  assumption  was  made  by  Ahrens  (1977)^  to  predict  irregular  wave 
runup  and  appears  reasonable  but  possibly  somewhat  conservative;  i.e., 
it  will  probably  yield  higher  extreme  vali  es  of  runup  than  are  actually 
observed.  The  procedure  for  combining  th  assumption  of  a Rayleigh 
distribution  of  wave  runup  with  equation  (1)  for  monochromatic  wave  over- 
topping is  outlined  below. 

1.  The  Rayleigh  runup  distribution  is  given  by 


( 


In  (1/p) 
2 


(2) 


^AHRENS,  J.P.,  "Prediction  of  Irregular  Wave  Runup,"  CETA  77-2, 
U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center, 
Fort  Belvoir,  Va. , July  1977. 
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where  Rp  is  the  wave  runup  associated  with  a particular  probability  of 
exceedance,  p;  Rg  is  the  wave  runup  with  a given  deepwater  significant 
height  and  period;  and  is  the  ratio  of  R-  to  Rg,*  the  index  "i"  is 

1 ^ — t t 1 , ___  _ 


employed  for  use  in  computations  shown  below.  Rg  is  referred  to  as  the 
significant  runup  and  is  predicted  in  the  conventional  manner  based  on 
laboratory  tests  using  monochromatic  waves  as  shown  in  the  SPM.  Equation 
(2)  provides  the  transition  between  monochromatic  and  irregular  conditions. 

2.  Equation  (1)  can  now  be  rewritten  as 


0.217 


h - d. 


where  is  the  overtopping  rate  associated  with  a particular  value  of 

Kj- , and  is  the  equivalent  deepwater  significant  wave  height. 

(H^)g  has  replaced  in  equation  f3)  to  emphasize  the  distinction  be- 

tween this  equation  and  the  monochromatic  overtopping  equation  (eq.  1), 
and  denotes  that  it  is  the  significant  height  which  is  being  used  in  this 
procedure.  The  term  (h  - d„)/R„  will  be  referred  to  as  the  relative  free- 
board. ® ^ 

3.  The  average  overtopping  rate  for  runups  with  a Rayleigh  distribu- 
tion, Q^,  was  computed  by  equation  (3)  using  the  following  relation: 


■4E 


p = 0.005  X i,  i = 1,  2,  3 . . . 199  . 

The  index  "i"  causes  a range  of  probabilities  to  be  generated;  these 
probabilities  when  substituted  into  equation  (2)  generate  a range  of 
values  of  the  parameter  K^.  K^,  in  turn,  is  substituted  into  equation 

(3)  to  generate  the  Q-^'s.  The  199  values  of  overtopping  are  computed 
for  a fixed  value  of  the  relative  freeboard  corresponding  to  the  various 
values  of  which  are  generated  by  varying  the  probability  of  exceed- 

ance, p,  from  0.005  to  0.995  in  increments  of  0.005.  The  average  value 
of  the  Qi's  (eq.  4)  represents  the  average  rate  of  overtopping  expected 
for  a time  interval  of  199  waves. 

4.  Using  the  procedure  shown  in  steps  2 and  3 above,  the  largest 
value  of  corresponds  to  a value  of  p = 0.005.  Since  the  peak  rates 

of  overtopping  might  be  of  concern  in  some  design  situations,  the  follow- 
ing definition  is  useful: 

Qq.5  = associated  with  a probability  of  exceedance 
of  p = 0.005  or  0.5  percent  . 
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If  equation  (2)  is  solved  for  p = Q.QQ5,  = 1.628  indicating  that 
Qq  g is  the  overtopping  rate  for  a runup  about  63  percent  greater  than 
Rg."  Qq  ^ will  be  referred  to  as  the  peak  overtopping  rate. 

The  Appendix  is  a FORTRAN  subroutine  which  calculates  Q^,  Q^,  5, 

Qp/Qf,.,  and  to  illustrate  how  the  above  procedure  can  he  program- 

ed and  applied  to  any  situation  which  might  occur. 

**************  example  problem  1************** 

This  example  is  based  on  the  example  worked  in  SPM  (Sec.  7.22), 
modified  to  account  for  the  irregularity  of  natural  wave  and  runup  con- 
dit ions . 


GIVEN : An  impermeable  structure  with  a smooth  slope  of  1 on  2.5  is 
subjected  to  waves  having  a deepwater  height  = 5 feet  and  a period 
T = 8 seconds.  The  depth  at  the  structure  toe  is  dg  = 10  feet  and  the 
structure  elevation  is  5 feet  above  Stillwater  level  (SWL) . Onshore 
winds  of  35  miles  per  hour  are  assumed. 


F IND : Estimate  the  overtopping  rate  for  the  given  wave. 


SOLUTION : The  wave  height  and  period  given  are  assumed  to  be  the  sig- 
nificant values;  therefore,  the  runup  obtained  in  the  example  is  the 
significant  runup,  R^.  Using  the  following  values,  as  given  in  the 
example  in  SPM, 


and 


a =0.08 


= 0.035 

H'  = (H')g  =5.0  feet 


solving  equations  (3)  and  (4)  for  Qj,  gives 

= 2.70  cubic  feet  per  second  per  foot 
as  compared  to  the  value  obtained  in  SPM  of 


(^  = 5.1  cubic  feet  per  second  per  foot 

for  monochromatic  waves.  Thus,  the  irregular  wave  overtopping  is 
about  one-half  of  the  value  given  by  the  monochromatic  wave  procedure. 
The  peak  value  of  the  irregular  wave  overtopping  is 


Qq  5 = 7.2  cubic  feet  per  second  per  foot  , 
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which  is  about  40  percent  greater  than  the  monochromatic  rate.  The 
wind  effect  on  overtopping  shown  in  the  SPM  example  is  not  discussed 
here,  but  it  can  be  assumed  that  it  would  increase  the  irregular  over- 
topping by  the  same  percentage  (11  percent)  as  the  monochromatic  over- 
topping. 

IV.  GENERALIZING  THE  IRREGULAR  WAVE  OVERTOPPING  TECHNIQUE 

Example  1 used  a specific  value  of  freeboard  for  the  structure 
(i.e.,  h - dg  = 5 feet);  however,  the  problem  could  have  been  worked 
for  any  value  of  freeboard  from  0 to  17  feet,  which  is  equivalent  to 
0.0  (h  - d3)/Rg  < 1.0.  The  example  can  he  considered  in  a more  gen- 
eral way  by  plotting  the  ratio  of  the  average  irregular  wave  overtopping 
rate,  Q^,  to  the  monochromatic  overtopping  rate,  (^,  as  a function  of 
the  relative  freeboard  (see  Fig.  1).  The  ratio  of  the  peak  to  monochro- 
matic overtopping  rates,  Qo.s/Qm*  also  shown  in  Figure  1.  When  the 
ratio  of  overtopping  rates  is  used  as  in  Figure  1,  the  influence  of  QJ 
cancels  out  and  the  curves  are  functions  only  of  a and  the  relative 
freeboard.  Figure  1 shows  the  curves  for  a = 0.08,  as  given  in  Example  1, 
and  by  reading  the  values  off  the  curves  for  a relative  freeboard  of 
0.294  the  results  of  the  example  can  be  verified. 

Figure  2 shows  a family  of  curves  of  Qp/%  and  Qq^^/Q^  for 
a = 0.04,  0.06,  0.08,  and  0.10  and  demonstrates  how  a controls  the  shape 
of  the  curves.  This  range  of  a's  includes  most  of  the  a's  tabulated  in 
the  SPM  (Sec.  7.22).  Figure  2 also  shows  that  generally  is  signifi- 
cantly larger  than  for  the  range  of  a's  presented.  The  irregular 

wave  overtopping  characteristics  shown  in  Figure  2 are  consistent  with 
the  trends  observed  by  Tsuruta  and  Goda  (1968)**  in  laboratory  tests  with 
both  monochromatic  and  irregular  waves. 

Figure  2 can  be  used  to  compute  for  various  values  of  a and 

(h  - dg)/Rg  by  regarding  the  ratio  Qy./C^  as  a "correction"  factor  for 
0^  as  illustrated  in  Example  2.  Since  it  is  relatively  difficult  to 
read  the  values  from  Figure  2,  the  ratios  of  Q^/Q^  and  Qg  5/Qm  are 
given  as  a function  of  a and  (h  - dg)/Rg  in  tabular  form  iii  Taoles  1 
and  2,  respectively. 

**************  example  PROBLEM  2************** 

To  illustrate  the  use  of  the  generalization  presented,  reconsider 
Example  1,  except  assume  the  freeboard  is  8 feet  rather  than  5 feet; 
therefore,  (h  - dg)/Rg  = 0.471.  Using  equation  (1)  for 

a =0.08 
05  = 0.035 

‘*TSURUTA,  S.,  and  GODA,  Y.,  "Expected  Discharge  of  Irregular  Wave 
Overtopping,"  Proceedings  of  the  11th  Conference  on  Coastal  Engineering, 
1968,  pp.  833-852. 


1 1 


nd  Qo.5  / Q 


0 2 ' ‘ I I J 

0 0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

Relofive  Freeboord  (h-ds/Rs) 


Figure  1.  Q^/Qm  Qo  5^^  curves  as  a function 

of  the  relative  freeboard  for  a = 0.08. 
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Relotiwe  Treeboord  (h-ds/Rs) 


Figure  2.  and  Qq.s/Qw  curves  as  a function  of  the 


relative  freeboard,  for  o » 0.04,  0.06,  0.08,  and  0.10. 


and 


% = (H^)s  = 5.0  feet 


0.471 


gives  a monochromatic  overtopping  rate  of 

0^  = 2.963  cubic  feet  per  second  per  foot  . 

Using  Table  1 for  a = 0.08  and  interpolating  between  (h  - dg)/Rg  = 0.45 
and  (h  - dg)/Rg  = 0.50  to  find  the  value  of  Qj,/C^  associated  with 
(h  - dg)/Rg  = 0.471  yield  a value  of  Qp/Qm  - 0.393;  therefore, 

= 0.393  X 2.963  = 1.16  cubic  feet  per  second  per  foot  . 

Using  Table  2 in  a similar  way  to  obtain  a peak  overtopping  rate  gives 

Qq  g = 1.784  X 2.963  = 5.29  cubic  feet  per  second  per  foot  . 

Thus,  Tables  1 and  2 give  factors  which  can  be  multiplied  by  the  mono- 
chromatic overtopping  rate,  given  by  equation  (1),  to  predict  the  average 
and  peak  overtopping  rates  for  irregular  wave  conditions.  In  addition, 
Tables  1 and  2,  along  with  Figure  2,  give  perspective  on  the  degree  of 
conservatism  or  nonconservatism  in  the  traditional  monochromatic  approach 
to  wave  overtopping. 

Example  problem  2 illustrates  the  output  format  of  the  computer 
program  given  in  the.  Appendix. 


V.  SUMMARY 

A proposed  technique  for  predicting  overtopping  rates  for  irregular 
wave  conditions  is  presented.  This  technique  is  easy  to  use  because  it 
can  be  applied  as  a "correction"  to  the  method  for  monochromatic  waves 
shown  in  the  SPM  (U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering 
Research  Center,  1975^;  Weggel,  1976^).  The  technique  has  the  further 
value  of  providing  insight  on  the  degree  of  conservatism  or  nonconserva- 
tism in  the  SPM  monochromatic  wave  overtopping  equation  (eq.  7-6,  eq.  1 
in  this  report)  for  different  wave  and  structure  conditions.  Insight  of 
this  type  can  be  used  in  evaluating  the  trade  offs  among  structure  height, 
structure  characteristics,  and  cost  when  determining  the  optimum  struc- 
ture design. 

Based  on  the  analysis  developed  in  this  report,  equation  (1)  generally 
predicts  overtopping  rates  substantially  greater  than  would  be  expected 


^U.S.  ARm',  CORPS  OF  ENGINEERS,  COASTAL  ENGINEERING  RESEARCH  CENTER, 
op.  cit. , p.  7. 

^WEGGEL,  op.  cit.,  p.  7. 
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If 


for  the  equivalent  irregular  wave  conditions,  except  for  structures  with 
high  relative  freeboards.  For  relative  freeboards  roughly  in  the  range 
0.85  < (h  - dg)/Rg  < l.Q,  the  irregular  wave  overtopping  rate  can  exceed 
the  monochromatic  rate.  In  addition,  there  will  be  a small  amount  of 
irregular  wave  overtopping  for  (h  - dg)/Rg  i 1-0  which  would  not  be  pre- 
dicted using  the  SPM  monochromatic  method.  The  analysis  also  shows  that 
the  peak  irregular  wave  overtopping  rate  (defined  as  Qq  will  exceed 
the  monochromatic  rate  for  all  values  of  the  relative  freeboard. 

Because  of  the  uncertainty  in  determining  the  wave  height,  runup, 
and  parameters  a and  required  to  predict  the  overtopping  rate,  the 

design  engineer  may  be  justified,  in  some  situations,  in  using  the  mono- 
chromatic wave  overtopping  technique,  which  is  generally  conservative. 

The  proposed  technique  should  yield  a more  accurate  estimate  of  the 
overtopping  rate  for  wind-generated  waves;  however,  it  is  not  necessarily 
conservative. 
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APPENDIX 

FORTRAN  Program  to  Compute 
Irregular  Wave  Overtopping  Rates 
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